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Abstract. The high sensitivity Parkes and Swinburne Multibeam pulsar surveys have nearly doubled the number of known 
pulsars and revealed many more distant pulsars with high dispersion measures. These data allow us to investigate in more 
detail the statistical parameters and distribution of pulsars, especially in the central regions of the Galaxy, which was almost 
impossible in previous low-frequency and less-sensitive surveys. To estimate the distances to pulsars we used the new NE2001 
Galactic electron density model. We study the Galactic distribution of normal pulsars with 1400 MHz luminosities greater than 
0. 1 mjy kpc 2 refining the shape and parameters of the radial distribution of pulsars. 

The maximum galactocentric distribution of pulsars is located at 3.2±0.4 kpc and the scale-length of this distribution is 3.8+0.4 
kpc for the assumed distances to the Galactic center R = 8.5kpc. The surface density of pulsars near the Galactic center is 
equal to or slightly higher than that in the solar neighborhood. 

For observable normal pulsars with luminosities > 0. 1 mjy kpc 2 , we also re-estimate their local surface density and birth-rate: 
41+5 pulsars kpc~ 2 and 4.1±0.5 pulsars kpc~ 2 Myr~' respectively. For the total number of potentially observable pulsars in 
the Galaxy, we obtain (24 ± 3) x 10 3 and (240 ± 30) x 10 3 before and after applying beaming correction according the Tauris 
& Manchester (1998) beaming model. Within the limits of errors of estimations these results are in close agreement with the 
results of the previous studies of Lyne et al. (1998) (hereafter LML98). 

The dependence of these results on the NE2001 model and recommendations for further improvement of electron density 
distribution are discussed. 
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Since the discovery of pulsars there have been numer- 
ous statistical investigations of their parameters (see for 
example Gunn & Ostriker 1970 Taylor & Manchester 1977 
|Lyneetal. 1985] IGuseinov & Yusifov 19841 |Narayan 
LML98, and references therein), such as the Galactocentric ra- 
dial and height distribution, the luminosity function, the evo- 
lution of magnetic field and luminosity of pulsars, etc. These 
results are then used for the estimation of the birth rate of pul- 
sars and for the modeling of Galactic structures. 

Early surveys of pulsars were carried out around fre- 
quencies of 400 MHz. At low frequencies the sky back- 
ground temperature dramatically increases in the direction of 
the Galactic center (GC), and also the interstellar medium 
(ISM) electron density rises. According to the sensitivity re- 
lation {Dewey et al. 1984*| see also relation (0} below) both 
of these effects lead to the reduction of search sensitivity to- 
ward the GC, which in turn, leads to the reduced detection 



Send offprint requests to: I. Yusifov, 
e-mail: yusifov@erciyes.edu.tr 



of pulsars around the GC. It makes it difficult to determine 
the precise radial distribution of pulsars in the Galaxy as well 
as the precise estimation of the number density of pulsars 
at the GC. So far in various statistical investigations, various 
approximations of radial variation of number density of pul- 
sars, such as filled center or hollow center distributions, were 
used (see for example |Narayan 1987 Hart man et al. 1 997 
IGonthier et al. 20021 [Schwarz & Seidel 2002 and references 
therein). A preliminary study of central regions of the Galaxy, 
on the base of early high frequency surveys (Cli fton et al. 1 992 
and lJohnston et al. 1992> . was carried out by Johnston (1994). 
The conclusion was the deficit of pulsars near the GC. 

Recent high frequency, sensitive Parkes Multibeam Pulsar 
Sur veys (PMPS) (Manchester et al. 200T1 IMorris et al. 20021 
and IKramer et al. 20031 revealed many more distant pulsars 
with high Dispersion Measures (DM). These data may al- 
low us to investigate in more detail the statistical parameters 
and distribution of pulsars, especially in the central regions 
of the Galaxy, which was almost impossible in previous low- 
frequency and less-sensitive surveys. Furthermore, high fre- 
quency searches reveal many pulsars at the lower as well as 
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X (kpc) 

Fig. 1. Projection of pulsars onto the Galactic plane. Pulsars discovered in the MBPS survey are marked by pluses; pulsars known 
before are marked by dots. GC is located at the beginning of coordinate system. Spiral arms of the Galaxy are from CL2002. The 
(x, y) coordinates of the Sun are (0.0, 8.5). The function of other lines and symbols is explained in the text. 



at the high luminosity ends of the luminosity function (LF) of 
pulsars ( Cami lo~2002l ILorimer 200*31 ). From the evolutionary 
standpoint it is of great interest in understanding the influence 
of such pulsars on the shape of pulsars' LF. 

We decided to estimate the density of pulsars at the GC and 
to take a fresh look at the radial distribution and luminosity 
function of pulsars using these new PMPS data from the ATNF 
Pulsar Catalogue of 1412 pulsars (Manches ter et al. 20 02 1. 
We estimated the distances to the pulsars using the NE2001 
Galactic electron density model, Cordes & Lazio (2002) and 
Cordes & Lazio (2003), (CL2002 and CL2003 hereafter). 

In the space and luminosity distribution studies, it is 
very important to make corrections for the observational se- 
lection effects. In various statistical studies, many authors 
(see for example [Lyne et al. 1985| |Narayan & Ostriker 1 990 
Bhattacha rya et al" 1992| ILorimer et al. 19931 IJohnston 19941 
Sturner & Dermer 19961 ICordes & Chernoif 19971 LML98, 
Gonthier et al. 200"2l and references therein) applied various re- 
liable but complicated statistical methods to precisely quantify 
important selection effects. In this paper we use an empirical 
(and simple) method for the correction of observational selec- 
tion effects, as in Kodiara (1974), Yusifov (1981), Leahy & Wu 
(1989) and Wu & Leahy (1989) (WL89, hereafter), with some 
modifications. 

A summary of the available data and apparent distribution 
of pulsars in the Galaxy are presented in Section 2. In Section 3 
we will discuss the selection effects and correction methods we 
applied. In Section 4 we derive the radial distribution of pulsars 



in the Galaxy. The discussion of the results and conclusions are 
presented in Sections 5 and 6. 

2. Available data and the apparent distribution of 
pulsars in the Galaxy 

At the time of preparing this paper, the number of exist- 
ing pulsars was 1412 (Ma nchester et al. 20021 ). The ATNF 
Pulsar Catalogue contains the data of PMPS, Swinburne 
( Edwards et al. 200T1 > survey and all previous pulsar survey re- 
sults, and it provides a good sample of data for the statistical 
study of pulsars. Nearly 600 of them were discovered in the 
PMPS survey. 

In this study we are interested mainly in statistics of "nor- 
mal" pulsars. For this reason we excluded from our sample bi- 
nary and recycled (P < 10~ 17 s/s), globular cluster and Large 
and Small Magellanic Cloud pulsars. Of course, some pulsars 
in the remaining sample are related to the recycled ones, but at 
the present stage it is difficult to separate them. This gives 1254 
pulsars in total, and 1043 of them are located in the regions of 
Galactic longitudes -100° < / < 50°. In this sample, 581 pul- 
sars were previously unknown and detected in the Parkes and 
Swinburne Multibeam Pulsar Surveys (MBPS). In our study 
we mainly used pulsars in the regions of the Galactic latitudes 
\b\ < 15° and longitudes -100° < / < 50°. However, intro- 
ducing the required corrections in density estimations, pulsars 
outside of this region were also taken into account, as described 
in the corresponding parts of the text. 
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Table 1. The apparent density of pulsars 



\R(kpc) 

\ 
r(Kpc) \ 





1.7 


3.4 


5.1 


6.8 


8.5 


10.2 


11.9 


13.6 


15.3 


17 


18.7 


0.0 












38.0 














1.7 










31.16 


16.12 


5.73 












3.4 








17.52 


15.14 


11.88 


4.06 


0.99 










5.1 






10.5 


11.63 


7.25 


6.75 


2.35 


0.6 


0.04 








6.8 




3.84 


8.22 


6.51 


5.87 


4.7 


1.69 


0.26 


0.07 









8.5 


1.29 


3.32 


3.04 


3.69 


3.41 


2.77 


2.13 


0.36 


0.1 


0.02 







10.2 




1.45 


1.69 


1.13 


0.8 


1.05 


1.27 


0.66 


0.09 


0.03 








11.9 






0.71 


0.51 


0.31 


0.32 


0.8 


0.42 


0.15 


0.09 


0.02 





13.6 








0.18 


0.11 


0.06 


0.23 


0.26 


0.2 


0.07 


0.07 


0.01 


15.3 










0.02 


0.13 


0.07 


0.22 


0.22 


0.18 


0.03 


0.01 


17.0 












0.06 


0.08 


0.04 


0.12 


0.2 


0.05 





18.7 














0.03 


0.01 


0.01 


0.04 


0.07 


0.01 



Projection of these pulsars on the Galactic plane 
is shown in Fig. ^ Although the recent observations 
( Eisenha uer et al. 20 03 1 show that the distance to the GC is 
8.0+0.4 kpc, in this study, as in the Galactic electron density 
model NE2001, we used the old value of the Sun-GC distance 
R Q = 8.5 kpc. 

The apparent distribution of pulsars for the subsequent cor- 
rections due to selection effects is derived in the following man- 
ner. We drew equidistant concentric circles around the Sun (r,) 
and the GC (Rj) and made up a quasi-regular grid of points at 
the points of intersection of these circles on the Galactic plane 
(see Fig.[0. To subdivide the Sun-GC distance into a sufficient 
number of intervals, distances between adjacent circles were 
selected equal to AR = R Q /5 = 1.7 kpc. Then we drew cir- 
cles of radius R C j around the grid points ( or cells ) and counted 
the number of pulsars within the boundaries of circle R C j. With 
the increase in distance from the Sun, the apparent density of 
pulsars decreases due to the selection effects. To avoid the loss 
of some pulsars at large distances from the Sun, the radius of 
the distant cells must be large. Arbitrarily, R C j increasing lin- 
early is chosen from Q.IR Q to 1/10 of the largest considered 
distances from the Sun (18.7 kpc) and calculated by the rela- 
tion R C j - 0.85 * (1.2 * (i - 1)/11 + 1), where i varies from 1 to 
12. Large cell radii at larger distances naturally reduce the den- 
sity fluctuations due to small number statistics and the selection 
effects. Grid circles around the Sun and GC were drawn until 
18.7 kpc. Apparent densities of pulsars (in units pulsars kpc -2 ) 
in grid cells are shown in Table 1 . 

Columns in Table 1 correspond to the cells at the same dis- 
tances Rj from the GC, which are indicated in the first row, and 
rows correspond to the cells at distances r, from the Sun, which 
are indicated in the first column. 

Assuming galactocentric symmetry of distributions, the av- 
erage densities in Table 1 are estimated in the following man- 
ner: for the regions |/| < 50° and |/| > 100°, the number of 
pulsars in symmetric cells around the y axis is averaged ( i.e. 
for the cells with coordinates ( y, x ) and ( y, -x ) in Fig.Q; for 
the region -100° < I < -50°, if the number of pulsars in the 
cells with positive x is more than in the negative x, we again 
averaged the numbers in the cells (y, x) and (y, -x). 



As is seen in Fig.^ at large distances from the Sun and the 
GC, the number of pulsars rapidly decreases, and in some cells 
pulsars are completely absent. In such cases, in order to reduce 
rapid fluctuations, in general, observational data are smoothed. 
Smoothing on the line of the Sun-GC direction, but behind the 
GC and in the anticenter direction, is carried out by the relation: 

Pij = (Pij- 1 + Pij+ 1 + 4 Pij + Pi-y + Pi+i/)/8, ( 1 ) 

where, p' t ■ is the smoothed surface density at the intersections 
of i* and / h circles correspondingly from the Sun and the GC. 
But in other directions the situation is changed, and each cell is 
surrounded by 6 other cells. For this reason smoothing for the 
distances more than 10 kpc from the Sun and GC (correspond- 
ing to the columns Rj > 1 1 kpc and rows r, > 11 kpc) is carried 
out by the relation: 

p\j = (Py-i +A-ij-i +Pi-ij + 6 Py+ 

Pi+lj +Pi/+1 + Pi+lj+1 )/12, (2) 

where designations are as in Eq. The densities in Table 1 
are obtained after applying these averaging and smoothing of 
observational densities. In these procedures, smoothed values 
are not additionally re-smoothed. For this reason, the cells to- 
tally surrounded by the empty cells in Table 1 remain with zero 
density. 

3. The selection effects and correction method 

Data from Table 1, after correcting for the selection effects, 
may be used for the calculation of the radial distribution of pul- 
sars in the Galaxy. 

In previous population studies (LML98) the number and 
birth-rate were estimated for the pulsars of 400 MHz luminosi- 
ties greater than 1 mJy kpc 2 . If we take into account the relation 
S ~ v" for pulsars, where S is observed flux, v is frequency 
and < a >= -1.7 ( Sieber 2002 1, the corresponding minimum 
luminosity at 1400 MHz is Lmi„i4oo ~0.1 mJy kpc 2 . Within a 
radius of 1 kpc from the Sun there are nearly 90 pulsars and 
only 10 of them have 1400 MHz luminosities < O.lmJy kpc 2 . 
In addition, for 4 of them flux densities at 1400 MHz are not 
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known and the approximate luminosities are estimated from the 
assumed mean spectral indices < a >= -1.7, so that the pop- 
ulation of low luminosity pulsars still remains uncertain due to 
small number statistics. For this reason, and for simplicity of 
comparison of similar parameters of early results we consider 
the sample of pulsars with 1400 MHz luminosities greater than 

£minl400- 

WL89 considered three types of selection effects, one of 
which is connected to the much more sensitive Arecibo sur- 
vey (Sto kes etal. 1 986 1 in comparison with other ones (Green 
Bank ( IStokes et al. 1986> . Jodrell Bank ( |Clifton & Lyne 1986) 
and Molonglo (Manch ester et al. 1978> V To get a homoge- 
neous sample, they excluded Arecibo pulsars in galactic lon- 
gitudes 42° < I < 60°. As we considered in general MBPS 
survey region of the Galaxy (-100° < I < 50°), we examine 
only two types of selection effects with some modifications. 

We divide selection effects into two categories and define 
them by the relation: 



p(r, R, l(r, R)) = K(l)K(r) Po (r, R, l(r, R)), 



where I is the Galactic longitude; r and R distances from the 
Sun and the GC; K(l) direction and K(r) distance-dependent 
selection factors; p„(r,R,l{r,R)) is the true and p(r,R,l{r,R)) 
is the observed distribution of surface densities of pulsars on 
the Galactic plane. Table 1 is a numeric representation of 
p(r,R,l(r,R)). 

K(l) is connected to the background radiation which leads 
to variations of the survey sensitivity with Galactic longitudes. 
Survey sensitivity is a function of several parameters, and in its 
simplest case may be described by the relation (see for example 
|Dewey et al.l984||Lyne et al. 1985) : 



-■fiS' 1 + 



sky 



DM \ 

1 + mJy, 

DM 



(4) 



where 5 m ; n (mJy) is the minimum detectable flux; S' (mJy) is 
the minimum detectable flux at the telescope beam center in the 
ideal situation; ft is a factor representing the reduction of sen- 
sitivity for the off-center scan of pulsars by a telescope beam; 
J S ky is sky background temperature; Tr is receiver noise tem- 
perature; DM is dispersion measure; and DMq is the charac- 
teristic dispersion measure depending on the antenna, receiver 
and other parameters of the survey. Sky temperature depends 
on observational frequencies and Galactic coordinates. The sky 
temperature map at 400 MHz has been determined by Haslam 
et al. (1982). Their results are analytically fitted by Narayan 
(1987) and have the form: 



' sky 400 



(/, b) = 25 + 



275 



[l + (//42) 2 ][l+(£/3) 2 ] 



K, 



(5) 



where I and b are the Galactic longitude and latitude, respec- 
tively. This sky temperature, then, may be scaled to the sky 
temperature of the observing frequency v (MHz) through the 
relation given by Johnston et al. (1992) in the form 



skyv 



/408MHz\ 2 - 6 

^sky400l I K - 



(6) 




-50 50 

Galactic Longitude, I, (deg) 



Fig. 2. Dependence of the correcting factor K(l) as a function 
of Galactic longitude. 

For the MBPS search frequencies we obtain: r s kyi374(/, b) = 
0.04 r s ky4oo(/, b). As in WL89, the direction-dependent correct- 
ing factor Kit) is defined by the relation: 



(3) K(l) = 1 + 



sky 1374 



4(/,* = 0) 



T R 



(7) 



where 7 , R =21°K is the noise temperature of the Parkes multi- 
beam receiver ( Manchester et al. 2001 ). Neglecting latitude de- 
pendence may lead to some overestimation of observed den- 
sities. However, taking into account that the majority of dis- 
tant pulsars from which radial distribution is derived are low- 
latitude pulsars, we may use Eq. @ as the first approxima- 
tion. The correction factor K(l) which was calculated from Eqs. 
(|51-0 is shown in Fig. 13 We assume that reduction of sensitiv- 
ity according to Eq. @, due to this selection effect, makes the 
observational number of pulsars in grid cells decrease. For this 
reason, the observational number (or density) must be corrected 
by this selection effect, dividing by K(l). In the anticenter di- 
rection and in the cell where the Sun located (r < 0.85 kpc), 
K(l) is assumed to be 1 . 

Kodiara (1974) gives an empirical method for correcting 
distance-dependent selection effects for the study of SNRs. 
Further, in Huang et al. (1980), Yusifov (1981), Guseinov & 
Yusifov (1984), Wu et al. (1984), Leahy & Wu (1989), WL89 
and Case & Bhattacharya (1989) this method was used to study 
the Galactic distributions of pulsars and SNRs. With some 
modifications, the same method is applied here to study the 
Galactic distribution of pulsars from MBPS survey results. 

In the flux-limited surveys for each distance (d) there is a 
minimum detectable luminosity defined by the relation 

^min = ^min^ > (8) 

where S m \ n is derived from Eq. and, for the PMPS survey, 
is ss 0.2 mJy. Pulsars with luminosities lower than L m ; n are un- 
detectable. Other reasons that reduce the sensitivity of pulsar 
surveys are the scattering and scintillation of pulsar radiation 
by the interstellar medium, pulse broadening, intrinsic varia- 
tion of pulsar radiation, off center scan of pulsars, etc. 

Pulse broadening, scattering and scintillation are strong 
functions of distances (see for example CL2003 and references 
therein). With increasing distances from the Sun each of these 
factors reduces the detected number (or density) of pulsars. 
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Table 2. Corrected densities of pulsars, in units pulsars kpc 



\R(kpc) 

\ 
r(Kpc) \ 


1.7 


3.4 


5.1 


6.8 


8.5 


10.2 


11.9 


13.6 


15.3 


17 


18.7 













41.5 














1.7 








80.62 


30.12 


10.6 












3.4 






83.86 


57.89 


41.54 


13.89 


3.39 










5.1 




93.05 


89.1 


50.54 


44.29 


14.87 


3.81 


0.24 








6.8 


63.01 


124.5 


90.66 


76.28 


57.94 


20.05 


3.11 


0.81 









8.5 


39.16 98.28 


84.9 


96.05 


83.5 


64.4 


47.74 


7.92 


2.19 


0.33 







10.2 


81.29 


89.65 


55.85 


37.45 


46.19 


53.7 


27.08 


3.48 


1.09 








11.9 




73.55 


48.92 


27.57 


26.79 


64.66 


32.65 


10.96 


6.48 


1.65 





13.6 






34.02 


19.91 


10.32 


35.91 


38.97 


28.94 


9.41 


10.26 


1.44 


15.3 








5.58 


42.4 


21.15 


62.58 


60.67 


46.45 


7.47 


2.66 


17 










36.9 


44.82 


24.57 


63.5 


102 


24.9 





18.7 












30.67 


11.41 


8.46 


39.34 


64.99 


6.77 



Here we make the simplifying assumption that the combined 
effect of these factors define K(r) in Eq. Q, which may be de- 
scribed by the exponential law as: 

K(r) = K x (r)K 2 (r) ■ • ■ K t (r) = e^e** ■ ■ ■ e~ c ' r = e~ cr , (9) 

where K\(r), ^(r) etc. are the distance-dependent correction 
factors relating to the pulse broadening, scattering and other 
selection effects. The quantitative estimation of each of these 
factors independently is difficult, but the combined effect of 
these factors leading to reduced detection of pulsars away from 
the Sun may be estimated empirically as described below. 

We assume that the surface density (SD) of pulsars is sym- 
metric around the GC, and considering a galactocentric circle 
with the radius R G , from Eq. Q we obtain: 

p(r,R Q )=K(l)K(r)p (r,R Q ), (10) 

where apparent densities p(r,R ) correspond to the R = 8.5 
kpc column in Table 1. According to the assumption p o (r,R ) 
is constant and K(l) is known. 









9., 












6 




6 - 




* , * 



2 4 6 8 10 12 14 16 18 

Heliocentric Distance, r (kpc) 



Fig. 3. Variation of the surface densities of pulsars on the galac- 
tocentric circle R = 8.5 kpc at various distances from the Sun 
(column 6 in Table 1). Apparent densities marked by circles 
and densities corrected due to the direction-dependent selec- 
tion effect (K{1) ) are marked by pluses. The dotted line is a 
LMS fitting of them 



The observed distribution of p(r, R Q ) against the heliocen- 
tric distance r is shown in Fig. [3] We assume that observed 
exponential decay of SD in Figgis connected to the distance- 
dependent correction factor K(r) in relations 0, © and J10I . 
Fitting the data in Fig.[3]with a simple exponent we can derive 
K(r) and p (r,R ). 

The approximate value of the local SD of pulsars p (r,R Q ) 
within r < 0.85 kpc may be estimated using all available data, 
taking into account the different sensitivity of surveys as a func- 
tion of the Galactic longitude and latitude. The method of es- 
timation is described in Yusifov & Kiiciik (2004). For the pre- 
liminary value of p a (r, R Q ) with 1400 MHz luminosities greater 
than L m ini4oo we obtain 38 pulsar kpc 2 marked by pluses in 
Fig. |3] A more precise value of p (r,R ) will be obtained ex- 
trapolating the fitting lines of distribution p(r, R ) as a function 
of r in Fig. [5] 

Fitting the corrected data in Fig.[3]by the least mean squares 
(LMS) method for the apparent distribution of SD, we can de- 
rive K(r) and p () (r, R Q ) (during the fitting, the point around 13.6 
kpc is not taken into account as an outlier): 

ln(p(r, R Q )/K(l)) = \n(p K(r)) = ln Po - cr, (11) 

where lnp = 3.73+0.16 and c = 0.362 + 0.017. The correction 
factor K{r) we define as 

K(r) = exv(-cr), (12) 

The apparent densities (py) in Table 1 are corrected by the 
relation using K(l) and K(r) from Eqs. Q and (I12> and 
corrected SD (p,y c ) of pulsars are given in Table 2. 

Extrapolating the fitted line in Fig. [3] for the local SD of 
pulsars with luminosities more than 0.1 mJy kpc 2 , we obtain 
po = 41 ± 5 pulsars kpc -2 . Although the obtained local SD of 
pulsars slightly exceeds the value derived by LML98 (30+6), 
within the errors of estimations the values are in good agree- 
ment. If we take into account some possible reasons leading to 
an increase of the apparent local SD of pulsars, discussed in 
section 5, the agreement is further strengthened. 
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Table 3. Errors of surface densities of pulsars estimated by the relation (12) 



\R(kpc) 

\ 
r(Kpc) \ 





1.7 


3.4 


5.1 


6.8 


8.5 


10.2 


11.9 


13.6 


15.3 


17 


18.7 















4.8 














1.7 










9.5 


4.7 


2.7 












3.4 








12.4 


9.0 


7.2 


3.9 


1.9 










5.1 






16.9 


15.6 


10.5 


9.5 


5.1 


2.5 


0.6 








6.8 




16.7 


25.1 


19.8 


17.3 


14.2 


7.5 


2.8 


1.4 









8.5 


15.9 


27.4 


24.5 


25.8 


23.1 


19.2 


15.8 


5.8 


3.0 


1.2 







10.2 




30.7 


31.9 


23.3 


18.0 


19.9 


21.4 


14.2 


4.8 


2.7 








11.9 






36.2 


27.7 


19.5 


18.7 


30.2 


20.1 


11.1 


8.4 


4.2 





13.6 








29.7 


21.5 


14.7 


27.7 


28.4 


23.8 


13.0 


13.5 


5.0 


15.3 










15.0 


40.7 


27.2 


47.5 


45.8 


38.9 


14.8 


8.7 


17 












50.6 


53.3 


37.5 


60.5 


77.5 


35.6 


0.0 


18.7 














59.1 


33.8 


28.0 


59.8 


76.6 


23.6 



4. Radial distribution of pulsars in the Galaxy 

By averaging the corrected SD of pulsars in Table 2 we can 
derive the SD of pulsars as a function of radial distance from 
the GC. 

As the selection effects vary with the direction and distance 
of the Sun, errors of corrected densities also vary from point to 
point. The errors of corrected densities in corresponding cells 
in general were estimated on the basis of Eq. Q- 

We assume that errors in calculating corrected densities p,- 7 - e 
are mainly connected to the errors in the correction factor K(r) 
and apparent density py = Nij/AS,, where Njj and AS ; are the 
number of pulsars in the cell and the area of the cell (/, j). K(r) 
is an exponential function (Eq. H2\ ) similar to y — a exp(bx), 
where the error due to uncertainty in x is estimated with the 
standard relation ( |Bevington 1969) : o~ y /y = bcr x . In evaluating 
errors for Njj we assumed that the number of pulsars in cells, 
in a rough approximation, follows Poisson statistics, that is the 
mean value equals the square of the standard deviation (fi = 
o 2 ). Depending on these assumptions, errors for py c can be 
calculated by the relation 

^ = UrJ +(rcrc) ' (13) 

where cr Pi c is the error of py c , and cr c = 0.017 is the standard 
deviation of c from Eq. (II It . In evaluating errors for the cir- 
cumsolar region (where r — in Eq. dl3» we took into account 
the error for po in Eq. dl It . Errors of densities for each cell 
of Table 2 estimated by this method are given in Table 3. From 
Tables 2 and 3 we can estimate the mean value and errors of SD 
of pulsars. As the errors for various data points are not equal, 
we calculated the weighted average of SD and corresponding 
errors for every column (or corresponding distances from GC) 
of Tables 2 and 3 by the following equations: 

p(R) = jr^- and o 2 (K) = \Y(o%)\ , (14) 

^r 1 ijc ^ i ' 



where the summation is carried out by index i on every column 
j. The radius Rj from GC corresponding to column j is given 
in the first rows of Tables 1-3. During the calculation of p(R) 
and cr(R), the densities marked in bold in Tables 2 and 3 are 
not taken into account, since they remain outside the region 
-100° < / < 50°. 

From Table 2 we see that the SD of pulsars at the GC is 
nearly 40 pulsar kpc 2 . This is the average density of pulsars 
within the cell of radius 1.3 kpc around the GC. In the consid- 
ered sample, within the galactocentric radius of 0.5 kpc pulsars 
are absent, however, from Table 3 one can appreciate the error 
of the density estimation in this region as ~16 pulsars kpc" 2 . 
But within the 0.5 kpc < R < 1 kpc from the GC there are 4 
pulsars. Applying correction factors K(l) and K(r) to the ap- 
parent density of pulsars around the galactocentric radius 0.75 
kpc we obtain 50 pulsars kpc" 2 . These values are shown in 
Figs.|4]and|5]and are used in fitting the radial distribution. 

All the results of obtained radial distribution of SD of pul- 
sars with corresponding error bars are plotted in Figs. |4] and 
|5] To simplify the comparison with other results, densities in 
Fig. |3 are normalized to the surface densities at the solar cir- 



100 - 




2 4 6 8 10 12 14 16 18 

Galactocentric Distance, R (kpc) 



Fig. 4. Radial distribution of surface densities of pulsars. 
Circles and error bars are derived from Eq. (I14> . Fitting line 
corresponds to Eq. ( fTSl 
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cle. Radial distribution of the SD in Fig.@]may be fitted by the 
frequently-used Gamma function written as: 



'"Ml; 



exp 



-b 



R-R 



Re 



(15) 



where R Q = 8.5 kpc is the Sun-GC distance. However, the 
relation (1151 implies that the surface density at R = is zero, 
which may be inconsistent with the reality. To obtain nonzero 
density at R — we include additional parameter Ri and used 
a shifted Gamma function, replacing R and R Q in Eq. (1151 by 
X — R +R\ and X Q — R G +R\ correspondingly. The best results 
of the fitting by the LMS method are: A = 37.6 + 1.9kpc" 2 , a = 
1.64 ± 0.11, b = 4.01 +0.24 and Ri = 0.55 + 0.10 kpc. 

The obtained radial distribution (Fig. 0} shows a clear de- 
crease of the SD of pulsars towards the GC, relative to the re- 
sults of early studies at low frequencies ( |Lyne et al. 19 85 1. The 
existence of pulsar deficit at the GC was already concluded by 
Johnston (1994), analyzing early high frequency pulsar surveys 
( IClifton et al. 1992l and Pohnston et al. 1992i . MBPS allows to 
more precise determination of the density distribution of pul- 
sars around the GC. As was shown earlier, at 0.75 kpc from the 
GC the SD of pulsars is around 50+16 pulsars kpc 2 and is not 
less than that in the circumsolar region. Applying the Tauris & 
Manchester (1998) (TM98 hereafter) beaming model, for this 
value we obtain (500+150) pulsars kpc -2 . 

At the present time there is a lot of information on radial 
distributions of various types of Galactic components, and it 
will be very interesting to compare radial distributions of pul- 
sars and SNR with the radial distributions of their progenitors, 
shown in Fig. EI 

It is believed that, in general, NS are formed during SN II 
explosions and that their progenitors are OB type Population I 
stars. However, from Fig.|5]it is seen that, although the maxi- 
mums of Population I objects coincide, the maximum of pul- 
sar distribution is shifted towards the GC by nearly 1.5 kpc. 
The maximum of pulsar and SNR distributions nearly coincide, 
but the radial scale length ( RSL ) of the pulsar distribution 
(~ 4 kpc) is nearly half that of the SNR distributions. 

The radial distribution of the birth location of NS is also 
striking. During population simulations of NS in the Galaxy, 
for the radial distributions of their birth location, the follow- 
ing relation derived by Paczynski (1990) (see for example 
ISturner & Dermer 1996l and lGonthier et al. 2002> is used: 



R 



R 



( ' R{R)=a iR^r v \-RT 

4.5 kpc and a« 



where R 



exp 



(16) 

1.068. The corresponding dis- 
tribution is shown in Fig. |5]by dots and marked as P. Although 
the maximum of this distribution coincides with the maximum 
of the SNR distribution, its shape considerably deviates from 
the distributions of Population I objects. 

If the progenitors of NS are the OB type Population I stars, 
then from the qualitative considerations of Fig. [5] it seems that 
the radial distribution of the birth location of NS must be lo- 
cated near the line shown by pluses in Fig. [5] and may be de- 
scribed by the relation: 



exp 



-b\ — 

Uo 




(17) 



0.5 1 1.5 2 

Galactocentric Distance, R/R (kpc) 



Fig. 5. Radial distributions of pulsars (squares) and other types 
of Population I objects: SNR distributions are from Case & 
Bhattacharya (1998); H2 column densities from Bronfman et 
al. (1988) and Wouterloot et al. (1990); HII regions are from 
Pladini et al. (2002), (HII regions represent the number of 
sources in the 0.5 kpc wide Galactocentric rings, but not surface 
densities); OB star formation regions data from Bronfman et 
al. (2000); the radial distributions of birth location of NS, from 
Paczynki (1990), (Eq. (I16» are marked by P; the expected ra- 
dial distributions of birth location of NS (Eq. dl7» are marked 
by pluses 



where a — 4, b — 6.8 and R Q = 8.5 kpc is the Sun - GC 
distance. The constant A must be chosen from the calibration 
constraints (particularly, in the current case A = 1050). The fi- 
nal solution of radial distributions of birth locations of NS may 
be possible after performing a detailed analysis as in Paczynski 
(1990), taking into account the new data on pulsar velocities, 
radial distributions of pulsars and comparing obtained results 
with the new radial distribution data of Population I objects. 

5. Discussion 

From the derived distributions and age of pulsars one can esti- 
mate the total number of pulsars in the whole Galaxy, birthrate 
of pulsars and other statistical parameters. New age estimations 
of pulsars requires population synthesis or careful analysis of 
the P — P diagram and a new study of z-distribution with more 
precise estimated distances of pulsars. However, from Kramer 
et al. (2003) and from our analysis it is obvious that distance 
estimations of pulsars require further improvement and addi- 
tional studies. For these reasons, we estimated the birthrate of 
pulsars by assuming t = 10 7 yr as the mean lifetime of active 
pulsars, leaving detailed study to the future. 

The total number of pulsars in the Galaxy may be esti- 
mated by integrating the radial distribution of SD from Fig.|4] 
Integrating the relation (I15i for the number of potentially 
observable pulsars with 1400 MHz luminosities greater than 
0.1 mJy kpc 2 , we obtain (24 + 3) x 10 3 and (240 + 30) x 



8 



Yusifov & Kiifiik: Radial distribution of pulsars 



10 3 before and after applying the mean value of the beam- 
ing factor (BF) from the TM98 model. For the correspond- 
ing birthrates we obtain nearly one pulsar for every 400+50 
and 40+5 yr in the Galaxy which is consistent with the 
LML98 results (one every 330 and 60 yr). Last (corrected) 
values are in close agreement with the recently derived rate 
of core collapsed supernovae (2 per century) in the Galaxy 
by lMcKee & Williams (1997)land (0 024-0.027)(g o /75) 2 yr 1 
by jvariTleTHjer^Thr^lcaur^l 994)| In this relation H„ is the 
Hubble constant. 

At the end of section 3 the local SD of pulsars was esti- 
mated as 41+5 pulsars kpc 2 . The corresponding local birth- 
rate of pulsars will be 4.1+0.5 pulsars kpc~ 2 Myr -1 . While the 
new values of SD and birth-rate are higher, within the errors of 
estimations they are consistent with the LML98 results (30+6 
and 2.8+1.7 correspondingly). 

MBPS surveys revealed new pulsars not only in the cen- 
tral and distant regions of the Galaxy, but also some local faint 
pulsars in the circumsolar region. It seems likely that this is 
one of the reasons for the slightly high local density of pul- 
sars obtained. Another reason for the apparent increase of the 
local SD probably is overestimated electron densities of ISM. 
This implies that the derived local SD is an upper limit and 
that the true value of the local SD of pulsars with luminosities 
greater than 0.1 mJy kpc 2 is located between 41 and 30 pul- 
sars kpc~ 2 , in even better agreement with the value derived by 
LML98 (30+6). 

Applying beaming correction according to the TM98 
model, for the local SD and birth-rate of pulsars we obtain 
(410 ± 50) pulsars kpc 2 and 41 ±5 pulsars kpc 2 Myr -1 , re- 
spectively. These values differ from the corresponding values 
of LML98 (156 pulsars kpc 2 and 10 pulsar kpc 2 Myr 1 ) 
by nearly 4 times. This is due to the BF used. We applied 
the TM98 BF (~10) which is nearly 2 times higher than the 
Biggs (1990)| BF, used in LML98. 

In evaluating the statistical parameters of pulsars, the re- 
liable estimation of distances is very important. Until re- 
cently the Taylor & Cordes (1993) (TC93 hereafter) model 
has been widely used for this purpose. The PMPS survey 
reveals much more distant pulsars with high DM and the 
following estimations show that the distances for these pul- 
sars calculated by the TC93 method are located far out- 
side of the Galaxy (Manchester et al. 200lllMorris et al. 2 002 
IKramer et al~2003l . 

The new electron density model developed by Cordes & 
Lazio (2002) (NE2001) is free of these disadvantages, and 
there are no pulsars outside the Galaxy according to their 
model. Another advantage of their model is its ability to be 
easily modified as further ISM data is collected (voids, clumps, 
HII regions). 

The mentioned electron density models consist of several 
components. One of them denoted as n\h\, the "thick disc" 
component and due to its large scale-height {h\ as 900pc) it 
occupies a significant volume of the Galaxy. In the NE2001 
model this value (0.033) is nearly two times higher than that of 
the TC93 model (0.0165). 

Analyzing the NE2001 model, Kramer et al. (2003) dis- 
cussed the decreasing trends of mean z-height with increasing 



distances of pulsars from the Sun. But in the TC93 model the 
mean z-height increases away from the Sun. The most probable 
reasons for a similar variation of z-height may be related to the 
above mentioned 2 times higher value of n\h\ in the NE2001 
model in comparison with the TC93 model. At the same time 
this also leads to some increase of local density of pulsars, as 
noted above. Some intermediate value of n\h\ possibly may 
improve the variation of the mean z-height away from the Sun 
and the local SD of pulsars. However, this demands the fine 
tuning of the electron density distribution model in the Galaxy, 
which is beyond the scope of this paper. 

In the future, in constructing new electron-density models, 
besides the improvements recommended by the authors of the 
NE2001, one must take into account the z-distribution of dis- 
tant pulsars, new HI absorption measurements of distant pul- 
sars and other statistical or integral parameters of the Galactic 
distributions of pulsars, derived in previous sections. Some im- 
provements also may be achieved taking into account warp and 
flaring structure of the Galaxy (Yusifov 2003 1. This may be 
possible after HI absorption measurements of a sufficient num- 
ber of high DM pulsars. 

The dimensions of galaxies, RSL of stellar distribution, 
and the location of star formation regions are important pa- 
rameters of galaxies. To estimate these parameters, in paral- 
lel with radio and optical data, researchers have applied in- 
frared (IR) and other measurements of various galactic compo- 
nents. For example, Drimmel & Spergel (2001), by analyzing 
far-infrared (FIR) and near-infrared (NIR) radiation data from 
the COBE/DIRBE instrument, constructed a three dimensional 
model of the Galaxy and found the value 0.287? G for the RSL 
of the stellar disc. In previous estimations this parameter has 
varied between 0.2 and 0.94 R Q (see Drim mel & Spergel 200 1| 
and references therein). 

OB stars are progenitors of NSs, and the radial distribution 
formation regions therein are studied in detail by Bronfman et 
al. (2000). The radial distribution of pulsars and OB star for- 
mation regions are compared in Fig. [5] and despite the fact 
that the maxima do not coincide, the agreement is reason- 
able. In general, the majority of star formation regions and 
supernova (SN) activity are located in spiral arms and in the 
molecular ring around the GC. Recent investigations show 
that OB star-forming regions and SN activity may also be lo- 
cated at the periphery of the Galaxy (see S til & Irwin 200 H and 
|Uyaniker & Kofhes 2002) . Distant pulsars with high DMs may 
play a crucial role in studying peripheral regions of the Galaxy. 
Once again this underlines the importance of reliable estima- 
tions of pulsar distances, especially for distant pulsars. 

6. Conclusions 

We have studied the population of normal pulsars with the lu- 
minosities Lhoo > 0. 1 mJy kpc 2 in the Galaxy, on the basis of 
the 1412 ATNF pulsar catalogue where the distances are cal- 
culated according to the new electron density model NE2001. 
Our results are as follows: 

- We have refined the radial distribution of surface density 
of normal pulsars by using the new distances. The maximum 
of radial distribution is located at 3.2 ± 0.4 kpc and ~1.5 kpc 
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nearer to the GC relative to the maximum of distributions of 
Population I objects. Although the maximum of distribution 
nearly coincides with the maximum of distributions of their 
progenitors (SNR), the RSL of pulsars is 3.8 + 0.4 kpc, i.e. 
nearly two times less than the SNR distribution. Integrating the 
radial distribution for the total number of normal pulsars we 
obtained (24 + 3) x 10 3 and (240 + 30) x 10 3 before and after 
applying beaming corrections according the TM98 model. 

- The local surface density of pulsars is ((41 ±5) and ((410+ 
50) pulsars kpc -2 ) before and after beaming corrections. 

- The surface density of pulsars around the GC region 
(R*0.7 kpc) is (50±16) pulsars kpc 2 and (500+150) pulsars 
kpc" 2 before and after applying beaming corrections. 

- The mean value of birth-rate of pulsars in the Galaxy is 
nearly one pulsar for every 400+50 and 40+5 yr before and 
after beaming corrections. 

- We recommend a new relation (Eq. dl7l i) for the ex- 
pected radial distribution of birth location of NSs, which is 
more closely related to the radial distribution of Population I 
objects. 

Acknowledgements. We would like to thank R.N. Manchester and 
the Parks Multibeam Pulsar Survey team for making the parameters 
of new pulsars available on the internet prior to formal publication 
as the ATNF catalogue of 1412 pulsars. We gratefully acknowledge 
the anonymous referee, for very detailed comments and suggestions, 
which led to significant improvements of the paper. We thank R. 
Wielebinski, J.L. Han and F.F. Ozeren for reading the manuscript and 
for useful discussions. We thank Victor B. Cohen for help in prepar- 
ing the manuscript. This work has been partially supported by Erciyes 
University R/D project No. 01-052-1, Turkey. Extensive use was 
made of both the Los Alamos preprint archive and the ADS system. 

Note added in proof: 

Recently, a paper of Lorimer, D.R., preprint| astro-ph/0308501 1 on the 
Galactic population of pulsars has been published. Although for study 
of pulsars they used different method, within the limits of the esti- 
mated errors, the obtained radial distribution and the total number of 
pulsars in the Galaxy (25±2)xl0 3 are nearly the same as in this paper. 

References 

Bevington, Philip R. 1969, Data Reduction and Error Analysis for the 

Physical Sciences, McGraw-Hill, New York, 63 
Bhattacharya, D., Wijers, R.A.M.J., Hartman, J.W. & Verbunt, E, 

1992, A&A, 254, 198 
Biggs, J.D., 1990, MNRAS, 245, 514 

Bronfman, L., Cohen, R.S., Ajvarez, H., May, J., Taddeus, P., 1988, 
ApJ, 324, 248 

Bronfman, L., Casassus, S., May, J., Nyman, L.-A, 2000, A&A, 358, 

521 

Camilo, F, 2002, preprint |astro-ph/02 10620 1 

Case, J.L. & Bhattacharya, D., 1998, ApJ, 504, 761 

Clifton, T.R. & Lyne, A.G., 1986, Nature, 320, 43 

Clifton, T.R., Lyne, A.G., Jones,A.W., McKenna, J., Ashworth, M., 

1992, MNRAS, 254, 177 

Cordes, J.M. & Lazio, T.J.W., 2002, preprint|astro-ph/0207156 1 

(CL2002) 

Cordes, J.M. & Lazio, T.J.W., 2003, preprint|astro-ph/0301598 1 

(CL2003) 

Cordes, J.M. & Chernoff, D.E, 1997, ApJ, 482, 971 



Dewey, R.J, Stokes, G.H., Segelstein, D.J., Taylor, J.H., Weisberg, 

J.M., 1984, in: Millisecond Pulsars, S.P Reynolds, D.R. 

Stinebring (eds.), NRAO, p.234 
Drimmel, R., & Spergel, D., 2001, ApJ, 556, 181 
Eisenhauer, F, Schodel, R., Genzel, R., Ott, T, Tecza, M., Abuter, R., 

Eckart, A. & Alexander, T, 2003, ApJ, 597, L121 
Edwards, R.T., Bailes, M., van Straten, W., & Britton, M.C., 2001, 

MNRAS, 326, 358 
Gonthier, PL., Ouellete, M.S., Berrier, J., O'Brien, S., & Harding, 

A.K., 2002, ApJ, 565, 482 
Gunn, J.E. & Ostriker, J.P, 1970, ApJ, 160, 979 
Guseinov, O.H. & Yusifov, I., 1984, Astron. Zh., 61, 708 (in Russian) 
Hartman, J.W. et al., 1997, A&A, 322, 477 
Haslam, C.G.T., Salter, C.J., et al„ 1982, Astr. Ap. Suppl., 47, 1 
Huang, K.L., Peng, Q.H., Huo, X.T. & Tong, Y, 1980, Acta Astron. 

Sin., 3, 237 

Johnston, S., Lyne, A.G., Manchester, R.N., Kniffen, DA., D'Amico, 

N, Lim, J., Ashworth, M., 1992, MNRAS, 255, 401 
Johnston, S., 1994, MNRAS, 268, 595 
Kodiara, K., 1974, PAS Japan, 26, 255 

Kramer, M., Bell, J.E, Manchester, R.N., et al., 2003, MNRAS, 342, 
1299 

Leahy , D.A. & Wu, X., 1989, PASP, 101, 607 

Lorimer, D.R., 2003, preprint|astro-ph/0301327 

Lorimer, D.R., Bailes, M., Dewey, R.J., & Harrison, PA., 1993, 

MNRAS, 263, 403 
Lyne, A.G. , Manchester, R.N. & Taylor, J.H., 1985, MNRAS, 213, 

613 

Lyne, A.G. , Manchester, R.N., Lorimer, D.R., Bailes, M., D'Amico, 
N, Tauris, T.M., Johnston, S., Bell, J.E & Nicastro, L., 1998, 
MNRAS, 295, 743 (LML98) 

Manchester, R.N., et al., 1978, MNRAS, 185, 409 

Manchester, R.N., Lyne, A.G., Camilo, F, et al., 2001, MNRAS, 328, 
17 

Manchester, R.N., et al., 2002, ATNF Pulsar Catalogue, 

http://www.atnf.csiro.au/research/pulsar/psrcat 
McKee, C.F & Williams, J.P., 1997, ApJ, 476, 144 
Morris, D.J., Hobbs, G., Lyne, A.G., et al., 2002, MNRAS, 335, 275 
Narayan, R., 1987, ApJ, 319,162 
Narayan, R. & Ostriker, J.P, 1990, ApJ, 352, 222 
Ord, S.M., Bailes, M. & van Staraten, W„ 2002, MNRAS, 337, 4090 
Paczynski, B., 1990, ApJ, 348, 485 

Paladini, R., Davies, R., DeZotti, G., 2002, 

preprint| astro-ph/02 1 234 1 
Schwarz, D.J. & Seidel, D., 2002, A&A, 388, 483 
Sieber,W., 2002, preprint|astro-ph/0208571 
Stil, J.M. & Irwin, J.A., 2001, ApJ, 328, 3591 

Stokes, G.H., Segelstein,D.J., Taylor, J.H., & Dewey, R.J., 1986, ApJ, 
311,694 

Sturner, S.J. & Dermer, CD., 1996, ApJ, 461, 872 
Tauris, T.M. & Manchester, R.N., 1998, MNRAS, 298, 625 (TM98) 
Taylor, T.M. & Cordes, J.M.,1993, ApJ, 411, 674 (TC93) 
Taylor, J.H. & Manchester, R.N., 1977, ApJ, 215,885 
Uyamker, B. & Kothes, R., 2002, ApJ, 574, 805 
van den Bergh, S. & McClure, R.D., 1994, ApJ, 425, 205 
Wouterloot, J.G.A., Brand, J., Burton, W.B. & Kwee, K.K., 1990, 
A&A, 230, 21 

Wu, X.J., Qiao, G.J., Dong, Z.W. & Lu, YY 1984, Pub. Beijing 

Astron. Obs., 7, 6 
Wu, X. & Leahy, D.A., 1989, Acta Ap. Sinica, 9, 233 (WL89) 
Yusifov, I., 1981, Former Soviet Astr. Circ, 1 164, 1 
Yusifov, I., 2003, in The Magnetized Interstellar Medium, eds. B. 

Uyamker, W. Reich & R. Wielebinski, Antalya, Turkey, pp. 165. 
Yusifov, I. & Kiiciik, I., 2004, in Preparation. 



